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,-. . Abstract 

qv \ We prove that the set of visible points of any lattice of dimension n > 2 

ON ' has pure point diffraction spectrum, and we determine the diffraction spectrum 

,^ ■ explicitly. This settles previous speculation on the exact nature of the diffrac- 

tion in this situation. Using similar methods we show the same result for the 
1-dimensional set of fcth-power-free integers with k > 2. Of special interest 
is the fact that neither of these sets is a Delone set -- each has holes of un- 
bounded inradius. We provide a careful formulation of the mathematical ideas 
underlying the study of diffraction from infinite point sets. 



Introduction 

It has long been known that the diffraction spectrum of a crystal consists of pure 
Bragg peaks only, being a pure point measure supported on the lattice dual to the 
lattice of periods. Until about 15 years ago it was tacitly believed that crystals were 
the only discrete point sets with this property. Now, however, we know that many 
quasicrystals have pure point diffraction spectra (though in a significantly different 
sense, since the locations of the peaks are no longer discrete) . These quasicrystals are 
all Meyer sets, that is, sets S that are both uniformly discrete and relatively dense 
and whose difference sets A = S — S also have these properties. (For a discussion of 
Meyer sets, see fl20|). There has consequently been a feeling that a perfectly diffracting 
discrete point set, if not precisely a Meyer set, must be closely related to a Meyer set. 
Indeed, the Meyer condition cannot possibly be strictly necessary since, as we shall see 
later, adding or removing a set of density zero does not alter the diffraction spectrum 
of a discrete point set, and one can clearly destroy the relative denseness of any 
uniformly discrete point set by removing a set of density zero. 

In this paper we give some simple examples of perfectly diffractive discrete point 
sets that deviate much further from the Meyer properties than this; in fact we consider 
sets which, for arbitrarily large D, have a lattice of holes of inner diameter at least D. 
Such a set cannot differ from a Meyer set (or from a model set, which is a particular 
case of a Meyer set) only by a set of density zero. The sets comprising our examples 
are well known in number theory: they are the sets of visible (or primitive) points of 
a lattice in any dimension n > 2, see [|1|, |nj] and the front cover of |1[ for a picture of 
the case n = 2, and the 1-dimensional sets consisting of the fcth-power-free integers 



for k > 2, see @, §6.6]. 

As well as giving a rigorous derivation of the diffraction spectra that depends on 
explicitly calculating the autocorrelation of the point set, we precede it by a shorter 
derivation of the pure point parts of the spectra only, via the Fourier transform of the 
point set. This takes a similar line to previous attempts and uses a general result of 
A. Hof. It provides a quick, elegant way of calculating the discrete parts of the spectra 
which gives the correct results but does not have a full mathematical justification at 
present. 

The rigorous approach is a response to the history of the problem for the visible 
lattice points, described in 0. In particular, there is a clear disagreement between 



earlier results in [29] and in |2T| regarding the nature of the diffraction. This was 
partially resolved in favour of |H| by a calculation of the point part of the diffraction 
spectrum in terms of Dirichlet series and its comparison with a real optical experiment 



0. The older numerical calculations in | 2"D[ , using the fast Fourier transform, suffer 
from an insufficient resolution and are misleading. In this article, we give the definitive 
description of the nature of the diffraction spectrum by making the previous formal 



calculation of the pure point part pT], £| rigorous and by proving that there is no 
continuous part. 

The paper is organized as follows. We start by describing the results from number 
theory and related areas we need, then describe some basic properties of the set of 
visible points. Next we discuss the background material of Fourier transforms and 
autocorrelations needed for diffraction spectra. This is first framed in the language 
of tempered distributions and then connected to measure theoretical considerations 
which are a basic part of the theory of diffraction. We then give the short intuitively 
motivated method of computing the pure point part of the diffraction spectrum, both 
for the visible points and the /cth-power-free integers. The following two sections are 
devoted to an explicit calculation of the autocorrelation of the set of visible points, 
followed by a rigorous derivation of the diffraction spectrum. The ensuing section 
treats the set of /cth-power-free numbers in the same way, followed by an extensive 
outlook and a summary. 

Strictly speaking, neither the measure theoretical picture nor the intuitive ap- 
proach to computing the diffraction are necessary for the logic of the paper, and they 
could be omitted by the reader whose primary concern is verifying the mathematical 
validity of the results. But we would like to stress that the additional information 
provided in this article is needed to locate our results in the wider context of mathe- 
matical diffraction theory. 

Tools from number theory 

Here we set out some results we need from number theory and related areas. 

Notation 

We use the notation (/, m) = gcd(Z, m) to denote the greatest common divisor of two 
integers I and m. For integers d and x, d \ x means that d is a divisor of x. Summation 
conditions like Li d \ x" are to be interpreted as meaning that d runs through positive 
divisors of x only (even when x is negative). The divisor function a(m) (defined for 
m G Z + ) counts the number of positive divisors of m, so a{m) = J2d\m 1- 

Also, we will frequently use the O-notation for error estimates. For example, we 
say that a function /(r) is 0(l/r) if there is a constant c such that |/(r)| is bounded 
by c/r for r > 0, see |I|, Sec. 3.2] for details. 



Lattices 

A lattice in M. n is a set T of the form 

T = Z6!©---©Z6 n , (1) 

where {b±, . . . ,b n } is a set of n linearly independent vectors called a basis of T. 
We define vol(T) to be the volume of a fundamental region of the lattice, e.g. of 
{t\bi + • • • + t n b n | < ti, . . . ,t n < 1}. Consequently, vol(T) can be calculated as 
vol(r) = | det(6i, . . . , 6 n )| which turns out to be independent of the basis chosen. 
Every lattice is uniformly discrete and relatively dense in M. n and is a subgroup of M. n 
under vector addition. 

Proposition 1 Let T be a lattice and a an arbitrary vector (not necessarily in T). 
Let N(R) be the number of points x inT + a with \x\ < R. Then there are constants 
c\ and C2, depending only on T, such that, for all R > 0, 

\N(R)vo\(T)-v n R n \ < Cl R n - l + c 2 , (2) 

where v n is the volume of the unit ball in W 1 , i.e. v n = 7r n / 2 /r(l + |). 

PROOF: The translates of the fundamental region of T by vectors in T + a tile M. n . 
Let V be the total volume of those translates that meet the open ball -Br(O) and v 
the total volume of those translates that lie entirely inside this ball. Then the volume 
v n R n of the ball and the number iV(i?) vol (T) are both bounded above by V and below 
by v. So their difference is at most V—v, which is the total volume of the translates 
of the fundamental region that meet the boundary of the ball. If D is the diameter^ 
of the fundamental region, then this volume is at most v n (R + D) n — v n (R — D) n < 
2 n v n DR n -\ when R > D, and at most v n (R + D) n < 2 n v n D n } when R < D. The 
second of these estimates is obvious, while the first (when R> D) follows from 



(R + D) n - (R - D) n = J2(}R n - m (D rn -(-D) 

m=0 ^ ' 

' n 



< 



2 £ Q R ^D (3) 

m odd x ' 

2DR n ~ x V ( n ) = 2 n DR n ~ 1 . 



m odd 



^^The diameter of a bounded set S C M. n is the supremum of all distances between points of S. 



This gives the result with c\ = 2 n v n D and C2 = 2 n v n D n . D 

If T is a lattice and r is a nonzero real number then rT is also a lattice (a basis 
of rT is {rb\, . . . ,rb n }). When r is an integer rT C T and rT is a sublattice (and a 
subgroup) of T, of index r n . For a nonzero integer m and two points a and b in T we 
write a = b (mod ml") to mean that a — b G ml\ With this notation we have the 
following Chinese Remainder Theorem for a lattice T. 

Proposition 2 Lei T 6e a lattice and a 1 , a 2 , . . . , a r G Y. If m 1; m 2 , . . . , m r G Z + 
are chosen so that (m-,m k ) = 1, for 1 < j < k < r, then there is a point a G T such 
that the solutions of the simultaneous congruences 

x = a x (mod m 1 T), . . . , x = a r (mod m r Y) (4) 

are precisely the points x ET with 

x = a (mod (m 1 • m 2 • . . . • m r )T) . (5) 



PROOF: This follows by applying the Chinese Remainder Theorem for integers ||14 , 
Thm. 2.7.1] to each coordinate with respect to a basis of T. □ 

We define the content of a nonzero lattice point a? in a lattice T by 

cont(cc) := max{ I \ x E IT} . (6) 

If x is expressed in terms of a basis of T, x = Y^ x jbj, then cont(:r) = gcd(xi, . . . , x n ) 
(which is therefore independent of the particular basis chosen). For consistency and 
convenience, we define cont(O) = oo. For m G Z + and x G T, we have 

cont (mcc) = m ■ cont(a;) . (7) 

It is clear from @ that for x G T \ {0} 

cont(aj) < J^L, (8) 

where \x\ is the Euclidean length of x and L(T) is the length of the shortest nonzero 
vector in T. 

The set V = V(T) of visible points of a lattice T (also known as the primitive 
points of T) is 

V := {xeT\ cont(a;) = 1} . (9) 

In terms of a lattice basis, V consists of all points whose coordinates have no common 
divisor. These are precisely the lattice points that are visible from the origin, in the 



sense that the line segment joining them to the origin contains no other lattice point. 
When n = 1, V consists of two points, equidistant from 0, but otherwise V is infinite 
and, indeed, contains more than half the lattice points, as we shall see. For m G Z + 
we have 

mV = {x G T | cont(aj) = m} . (10) 

Density 

Let 5* be a uniformly discrete set of points in M n (i.e. there is a c > so that x,y G S 
with x 7^ y implies \x — y\ > c) . We say S has natural density T> = dens(S') if 

^ ^ ,ns \{x e S \ \x\ < R}\ m . . 

V R = V R (S) := - — 11 — > x> (11) 

H Hy ! v n R n V ; 

as R —> oo. The expression \T> — T> R \ (a function of R) is called the error terra for the 
natural density. For example, by Proposition [TJ, D = dens(r) exists for every lattice 
T and is equal to l/vol(r) (with error term 0(1/ R)). 

If dens(S") exists and if T is an orthogonal transformation then dens(T(S')) = 
dens(S'), but it is a failing of the natural density that it is not always true that 
dens(T(S")) = dens(S')/| det(T)| when T is a general affine transformation, see Ap- 
pendix. 

When we use the word "density" from now on we shall mean natural density. 

Power-free numbers 

For an integer exponent k > 1 the set F = F^ of kth-power-free integers is 

F := {n G Z | n is not divisible by d k for any integer d > 1} . (12) 

An equivalent characterization of the numbers n G F is that in their prime power 
factorization n = ip" 1 ^ 2 ' ' 'P a r r every exponent a,j is less than k. In close analogy 
with V, the case k — 1 is trivial, F\ consisting of just the two numbers ±1, but for 
k > 2, Fk is infinite and contains more than half the integers. These numbers have 



been studied for a long time, see p2| , |T9| for two sources relevant in our context. 



Inclusion-exclusion 

The Mobius function fi(m) is defined for m G Z + by 

1, when m = l, 

fj,(m) := ^ (— l) r , when m is a product of r distinct primes, (13) 

0, when m is divisible by the square of a prime. 

6 



It is multiplicative in the sense that fi(lm) = fi(l)fi(m) when (l,m) = 1 (and clearly 
fj,(lm) = when (l,m) 7^ 1). 

A number of inversion formula? and variants of the inclusion-exclusion principle 
can be expressed in terms of this function. The result we need in this paper is 

5>o = (i' l m= \' (14) 

^-^ 1 0, if m > 1, v ' 

d\m 



see 



14| , Thm. 6.3.1]. For a point x in a lattice V this shows that 



Xv( x ) : = Yl M ( 15 ) 

d|cont(a;) 

is the characteristic function of the visible points V of T (except that it is undefined 
when x = 0). 

Similarly, since d k \ x if and only if d divides the largest integer whose fcth-power 
divides x, 

Xf {x) := Y,M ( 16 ) 

d k \x 

is the characteristic function of the A;th-power-free integers (except for being undefined 
when x = 0). 

Dirichlet series 



The well-known Riemann zeta-function (see, for example, |l], Ch. 12] or [14, §9.2]) is 
defined for complex numbers s with Re(s) > 1 by 



cw==E^ = nfi-ir. < 17 > 



m=l 



-1 

m s J " J " V P s 



where the sum is called a Dirichlet series and the product (in which p runs through 
all positive prime numbers) is called an Euler product. The sum and product are ab- 
solutely convergent for s in the half-plane Re(s) > 1, but ((s) can be meromorphically 
continued to the whole of C, for example £(0) = —1/2. The only singularity of ((s) 
is a simple pole at s = 1 with residue 1. Using the Euler product, we see that the 
Dirichlet series of 1/C( S ) is given by 

1 -nfi-=) = £* («) 



C(s) *■* V p s J ^ m 



m=l 



This function has infinitely many poles (all in the half-plane Re(s) < 1). Its value at 
1 is and its value at is —2. 

Another Dirichlet series we shall encounter is 

m=l P 

Again, both the Dirichlet series and Euler product are absolutely convergent in the 
half-plane Re(s) > 1. It can also be seen from the Euler product that £(1) = 0. 

Visible points of a lattice 

Here we summarize some elementary and well-known properties of the set of visible 
points V of a lattice T, together with complete proofs. 

Since T is a free Abelian group of rank n, its automorphism group, Aut(T), is 
isomorphic to the matrix group GL(n, Z). Explicit isomorphisms can be found by 
taking coordinates with respect to any lattice basis. 

Proposition 3 The orbits of the action of GL(n, Z) on T are the sets mV , m G Mo. 
In particular, GL(n,Z) acts transitively on V. 

PROOF: Since the elements of GL(n, Z) cannot decrease content and are invertible, 
they preserve content. Hence each set mV is invariant under GL(n, Z). The transi- 
tivity of GL (n, Z) on V can be seen from the facts that every visible point belongs to 
some basis of T 0, §3,Thm. 5] and that any two bases of T are related by a trans- 
formation in GL(n, Z). For m £ Z + the transitivity of GL(n, Z) on mV follows from 
its transitivity on V. Transitivity on the singleton orbit • V — {0} is trivial. □ 

Proposition 4 V is uniformly discrete, but has arbitrarily large holes. Moreover, 
for any r > 0, there is a set of holes in V of inradius at least r whose centres have 
positive density. 

PROOF: The uniform discreteness is trivial, as V is a subset of a lattice. Now let 
C = {a l5 . . . , a s } be any finite configuration of points in T (e.g. all points in a ball or 
a cube). Choose s integer moduli mi, . . . ,m s that are > 1 and coprime in pairs (for 
example, they could be the first s primes). By Proposition |2| there is a point a e T 
with 

a = — a x (mod m-^T) , . . . , a = — a s (mod m s T) . (20) 



Now for any x = a (mod vfi\Tfii . . . m s T) the configuration C+x = {a 1 +x, . . . , a s +x} 
is congruent, in the geometric sense, to C but no point in C + x is visible, since 
a,j + x G rrijT. The points a; have density dens(r)/(m 1 m 2 . . . m s ) > 0. D 

For n — 2, and with the density of holes not mentioned, this is Thm. 5.29 of [l|]. 
We note that the hole nearest the origin provided by this argument can be expected to 
be at a distance of the order s s and that the density guaranteed for holes of inradius 

2 n 

r is of the order r~ n r , so large holes, while having positive density, are probably 
extremely sparse. 

This proposition shows that V, though uniformly discrete, is not relatively dense, 
and hence not a Delone set. Consequently it is not a Meyer set either. (Recall that 
A is a Meyer set if and only if both A and A — A := {x — y \ x, y G A} are Delone 
sets [fL7| .) Also, V cannot be transformed into a Delone set by adding a set of zero 



density. However, we do have: 
Proposition 5 Ifn>2 then V — V = T. 

Proof: Clearly, V - V C T. Now let x = YTj=x x j b j e r > where n > 2. Then, 
(xi, %i-> ■ ■ ■ , x n ) = (xi + 1, 1, X3, . . . , x n ) — (1, 1 — X2, 0, . . . , 0), and x is the difference 
of two visible points. □ 

The existence of arbitrarily large holes also implies that the set of visible points 
cannot have a uniform density (not, at least, when its density is positive in some 
sense, as is the case for n > 2). It does have a natural density, however: 

Proposition 6 The visible points V of a lattice T G M n have a natural density given 
by 

d ™w = iSp- < 21 > 

with error term 0(1/R) when n^2 and 0(logR/R) when n = 2. 

This is a standard example of the use of Mobius inversion given (at least for the case 
r = Z 2 ) in most introductory number theory books (for example, [[H], Thm. 6.6.3] 
and [Q, Thm. 3.9]). In these particular references the averages are taken over triangles 
and squares, respectively, instead of balls. Indeed the density is independent of the 
shape of the region averaged over. We say more about this in the Appendix. 

PROOF: The proposition is trivially true for n — 1, when the pole of ((s) at 1 
gives a density of 0, so we assume from now on that n > 2. 

The density of V, if it exists, is the limit as R —>■ oo of 



^E 1 ' (22) 



v n R 

|a;|<_R 

xev 



which by ( [Op is 

^ E E*™> = ^ £ *E'. < 23 > 

ccer\{0} m|cont(as) l<m<cfl ccGmr\{0} 

|a=|<_R. |cc|<R 

where c = 1/L(T) (using (j8|)). The inner sum is equal to the number of nonzero points 
y ET with \y\ < R/m, which by Proposition [1] is 

§y+o((TfVod). (24) 



vol(r) \mj \ \m 



Substituting this into the right hand side of fl2"3| ) gives 



^ S f *"(i E ^t)^(^) (25) 

l<m<aR l<m<C-R 

which tends to dens(r)/C(n) as R — > oo when n > 2 by ( jI5|) . The total error is 
0(l/i? n_1 ), from the last term and the tail of the sum in the main term, and 0(1/ R) 
(or 0(logR/R) when n = 2 and the series diverges logarithmically), from the middle 
term. □ 

Calculations of densities by Mobius inversion form the core of this paper. This is 
the first of many. 

Diffraction spectra 

In this section we assemble the facts we need about distributions, Fourier transforms 
and diffraction spectra. The mathematics underlying diffraction is quite subtle and 
needs to be spelt out carefully. Although the discussion in this section does not contain 
much that is new, it is nonetheless difficult to extract it all from any convenient source. 
We shall use a formulation based upon tempered distributions. For an essentially 



parallel approach which starts with measures, we refer to ||12| . We shall link the two 
approaches later in this section. 

Autocorrelations 

In dealing with diffraction, and therefore Fourier transforms, it is appropriate to 
use tempered distributions, whose test space S consists of the Schwartz functions 
(also known as "rapidly decreasing functions"). We refer to ([|3D], 25 1) for this and 

10 



the details on the standard topology which is used to describe convergence in S. A 
tempered distribution T is a continuous (in the sense of this topology) linear functional 
on S. The space of tempered distributions is denoted by S' and is equipped with the 
weak*-topology. A simple example is S x , Dirac's delta- distribution^ at the point x, 
defined by (5 X , iji) := if){x) for all if; 6 S. We prefer this notation to 5 x {jp) because it 
emphasizes the duality between S' and S. 

If S is a uniformly discrete subset of R n , we call 



u s 



J2 6 * ( 26 ) 



xes 



its Dirac comb. It is also a tempered distribution, and in fact the sum (^) is conver- 
gent (with any ordering) in the weak*-topology. 

To describe the diffraction from a Dirac comb uo = u s we need its natural auto- 
correlation distribution (also called its generalized Patterson function) defined by 

7 - := &^ £ 5 *-^ (27) 

x,yeS R 

where Sr = S fl B R (0) and -Br(O) is the ball of radius R centre 0. We shall simply 
call 7^ the "autocorrelation" of S from now on. The existence of this limit in the 
weak*-topology is a prerequisite for the diffraction spectrum to be well defined. (The 
word "natural" refers to the use of the expanding ball B R (0) in the averaging process. 
Replacing it by an expanding region of some other shape might lead to a different 
limit.) It is clear from the definition that enlarging or diminishing 5* by a set of 
density does not change its autocorrelation j u . In particular, adding or removing 
any finite number of points does not change 7 W , as for natural density. 

We say that S has finite local complexity if A = S — S is discrete and closed in 
W 1 (as is clearly the case for all sets whose diffraction spectra we seek in this paper). 
Such an S is uniformly discrete because G A is isolated. The existence of the limit 
(p?1) in the weak*-topology is now controlled by the following result. 

Lemma 1 Let S be a set of finite local complexity and uo = us its Dirac comb. Then 
S has a natural autocorrelation 7^ if and only if the coefficients w(a) 

w(a) := lim V 1 (28) 

[sb|,|sc— a\<R 
x,x—a£S 



2 Nowadays it is usually called Dirac's point measure for reasons that will become clear shortly. 

11 



exist, i.e. the right hand side is convergent for all a G A. In this case, w(a) > for 
all a G A and 7^ is the tempered distribution of positive typ^ given by the weak*- 
convergent sum 

lu = £>(a)<J a . (29) 

aeA 

Proof: The existence of 7^ clearly implies the existence of w(a) for all a G A because 
A is discrete by assumption and the Schwartz space S contains all C^-functions of 
compact support, whence we can focus on any individual a. 

Conversely, assume that the w(a) exist (they are then clearly > 0). Since A 
is closed and discrete, its intersection with any compact subset of M n contains only 
finitely many points. Consequently, the w(a) are locally summable and the right hand 
side of (p9|) defines a distribution over the space V of all C^-functions of compact 
support. We have to show that it is actually also a tempered distribution. This follows 
from the translation boundednessQ of u which is then inherited by j u , see [|12|, Prop. 
2.2]. Finally, 7^ can also be written as a certain volume-normalized convolution (see 
below) which implies that it is a distribution of positive type. □ 

Fourier transforms 

The Fourier transform T of a tempered distribution T G S' is defined by (T,ip) : = 
(T, ip) , where we use the definition 

4>(y) := I e~ 2lTiy - x i)(x)dx (30) 

for the Fourier transform of functions ip G S. The Fourier transform maps the space S 
onto itself and is continuous on S in the standard topology for Schwartz functions p5 , 



Thm. 7.7], hence it maps S' onto itself and is continuous on S' in the weak*-topology 
2|, Thm. 7.15]. 

As special cases, we mention <5o = 1 and the well-known Poisson summation 
formula for lattice Dirac combs 

u T = dens(r) • u r , (31) 

where T* is the dual or reciprocal lattice defined by 

P := {y I y ■ x G Z , for all x G T} . (32) 

3 See next section for a definition. 

4 We will explain this below in the context of measures. 

12 



(Eq. pi| ) can easily be derived from Poisson's summation formula for Schwartz func- 
tions 0, p. 254].) 

For a set S of finite local complexity with autocorrelation j u , its diffraction pattern 
(also called its diffraction distribution or diffraction spectrum) is the Fourier transform 
j u . In view of the remarks above, 7^ is a tempered distribution. It is also a positive 
measure, as we shall see. 

The autocorrelation of a lattice T, for example, is supported on T itself (since 
A = T — T = T) and each peak has equal amplitude dens(r). So the autocorrelation 
is 7^ = dens(r)u; r and the corresponding diffraction spectrum is 7^ = dens(r) 2 co> r *, 
a constant multiple of the Fourier transform of ui T itself. However, as we shall see, 
the Fourier transform of a general point set (even of finite local complexity) does not 
describe its diffraction in such a simple way. 

Pure point distributions 

Our principal concern is with showing that the visible points and fci/i-power-free points 
have a pure point spectrum. In this section we consider a special class of point 
measures which we will use in the sequel. We already borrow from the terminology 
of measures here, although we will establish the precise connection only in the next 
section. 

Consider an arbitrary complex point measure. It can be expressed in the form 



^w(x)8 x , (33) 



xes 

where the point set S is countable, but not necessarily uniformly discrete, and the 
coefficients or weights w(x) e C are not necessarily constant. Note that the weights 
may be complex numbers. Let us assume in addition that the measure is translation 
bounded. This can be expressed as the condition that for every compact set K C M. n 
the sum 

£ W*)\ (34) 

xeSn(K+a) 

is convergent and bounded uniformly in a. We denote the space of all translation 
bounded point measures by T. All these measures are tempered, and we identify v 
with T„, the corresponding tempered distribution. 

As a subset of S', T is not closed in the weak*-topology. For example, the sequence 
of pure point distributions 

0'""<W ( 35 ) 
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tends to the constant function 1 as j —* oo. (For a test function ip the numbers 
j~ n uJ Zn ,Aip) are approximating sums to the integral^ of ip.) A similar argument shows 
that every bounded continuous function is a limit of pure point distributions, compare 



the more detailed discussion in |24], Sec. IV. 5]. Unfortunately, as this example shows, 
the weak*-limit of such a sequence is not, in general, the same as its pointwise limit. 
However, taking pointwise limits in T is justifiable under certain circumstances, 
and fortunately these apply in the cases of interest to us here. We introduce a "locally 
defined" norm on T by 



Mlioc : = sup / d\v\ = sup V |io(a?)|, (36) 

k Jk k jt£ v 



xesnK 

where the supremum is taken over all compact sets K of diameter < 1. This norm 
defines a topology on T stronger than the weak*-topology and it provides a simulta- 
neous "M-test" for pointwise and weak*-convergence of infinite sums of translation 
bounded point measures: 

Lemma 2 If Vj G T for j G Z + and YlT=i 1 1 ^i I hoc i> s convergent then YlTLi^j ^ s 
pointwise and weak* -convergent to the same sum v aT . 

Proof: Let Vj = J2 x &s w j(. x )$x and choose a fixed covering of {K m } meZ + of M. n 
by compact sets of diameter < 1. For any ip G S we have 

oo 

io^voi = | J2 w A x )i j i x ) < Yl J2 \ w j( x )^( x )\ 

x&Sj m=l x£SjC\K m 

oo 

< Ikjlhoc^ ll^xjl = C i>- INh°o (37) 

m=l 

where ip Km is the restriction of ip to K m and C^ < oo is a constant that depends 
only on ip and the covering chosen. Hence Yl^jLii^j^) * s absolutely convergent (by 
comparison with ^ ||^||i c) and so Yl u j i s weak*-convergent to a distribution v. 
Also 

oo oo oo 

(i/,v) = ^2(yj,ij>) = J2J2wj(x)ip(x) = Yl ^2 w 3( x M x ) > ( 38 ) 

i=i j=i xeSj xe\J Sj j=i 

the reversal of the order of summation (with Wj(x) := whenever x (jL Sj) being 
justified by the fact that the double sum is absolutely convergent, in view of (|37|) . 

5 Seen as a sequence of measures, ([55]) weak*-converges to Lebesgue measure. 
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Hence 

oo 

is the pointwise sum of the u/s. 

Finally, if K C W 1 is any compact set with diameter < 1 and we write S :— (J Sj 
and w(x) := Y^jLi w j( x ) then 



oo oo 



tceSn(A'+a) j=l xeSjn(K+a) 3=1 

Hence v is translation bounded, and therefore tempered. □ 

Note that if /ij is the distribution in (|35|) then fij +1 — \ij is translation bounded 
but 

Ik'+i - Millioc > ^TT - j(j + i) ' ( 41 ) 

where f n is the volume of the n-dimensional unit ball. (The supports of \Xj and Hj+i 
intersect in the integer points only and there is almost no cancellation in calculating 
the norm.) So the norms of the differences are bounded below by a positive constant 
for large j and their sum diverges to infinity, thus failing to satisfy the hypothesis of 
the lemma. (An easier calculation along the same lines is ||/i2 fc + 1 — A i 2 fc ll ^ v n /2 n .) 

Distributions and measures 

The diffraction pattern 7^ is a tempered distribution. However, it is also a positive 
measure. This remarkable fact is indispensible for the general theory of diffraction, 
making available to it a vast array of concepts and tools. For example, it makes it 
immediately evident that the diffraction pattern may be viewed as having a pure point 
part and a continuous part, a fact of considerable physical significance. 

The next two subsections sketch out this distribution theory - measure theory con- 
nection in the context of diffraction. The proofs of the main theorems of this paper in 
no way depend on this background material, which can therefore be ignored as far as 
verifying the results is concerned. Nevertheless, it is desirable to realize the natural 
connection to measure theory, both to see our results in their proper setting and be- 
cause it is this picture that has a generalization to the diffraction theory of translation 
bounded measures on locally compact Abelian groups which is an appropriate setting 
for more general questions, compare [IT, £|. 



A distribution T is positive if for all positive test functions ip, (T,ip) > 0. A 
distribution T is of positive type or is positive definite if for all test functions ip, 
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(T, ifj*ip*) > 0, where ip*(x) := ip(—x). These two concepts are related by the Bochner- 
Schwarz theorem (see |K], Thm. VII. XVIII] or [^4], Thm. IX. 10]) which asserts that a 
tempered distribution is positive if and only if it is the Fourier transform of a tempered 
distribution of positive type. 

Let C denote the space of complex- valued continuous functions of compact support 
on R n and let || ■ || denote the supremum norm on C. A (complex) measure v on M. n 
is defined as a linear functional on C such that for every compact set K C M. n there 
is a constant a K such that 

K0)l < a K U\\ (42) 

for all G C with support in K. Such measures are in one-to-one correspondence 
with regular Borel measures through the Riesz- Markov representation theorem, see 
|8|, Ch. XIII] and [|, Ch. 8, Sec. 69] for background material, and we thus identify 
these two pictures. We deal only with regular Borel measures in this paper. If for 
each <p G C inequality ( ^2] ) holds uniformly for all translates of <f), we say that v is 
translation bounded. This turns out to be a very useful concept because the Fourier 
transform of a tempered measure, though a tempered distribution, need not be a 
measure, but the Fourier transform of the autocorrelation of a translation bounded 
tempered distribution, if it exists, is always a measure. Let us explain how to control 
this subtle point in the context of diffraction theory. 

Distributions and measures act on different spaces of functions and are equipped 
with different topologies. Nonetheless, there is an important connection between them 
which comes through the fact that the space T> of C°°-functions of compact support 
is dense both in C and in S. 

If a measure v defines a tempered distribution T v by 

(T„,V) = v(i>) = fipdv (43) 

for all ip G S, the measure v is called tempered. A sufficient condition for a measure to 
be tempered is that it is slowly increasing in the sense that J(l + |a?|)~ fc |z/|((ia3) < oo 



for some k G Z + , see [30, Thm. VII. VII] or [25, Ex. 7.12 b]. Here, \u\ is the unique 
absolute value of v, i.e. the smallest positive measure p such that |^(0)| < p{\<t>\) for all 
G C. It is also called the variational measure of za As a partial converse, any positive 
tempered distribution is a positive tempered measure. Thus, under the assumption 
of positiveness, tempered measures and tempered distributions can be viewed as the 
same thing, and the Bochner-Schwartz theorem can be restated as follows: v is a 
positive tempered measure if and only if it is the Fourier transform of a tempered 
distribution of positive type. 
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Now, let v be a translation bounded measure. Clearly v is tempered. Our pre- 
vious definition of the autocorrelation of a Dirac comb has a natural extension to 
the autocorrelation r j v of v by means of a volume-normalized convolution of v with 
v*, where v* is defined by v*{<fi) '■— v{4>*), see [12[ for details^. This construction of 



7j, guarantees that, when it exists, it is a positive definite tempered measure (and 
distribution), so combining the previous arguments we see that the corresponding 
diffraction pattern j v is a tempered positive measure. 

This is important because it is this object that describes what one actually sees 
on the screen in a diffraction experiment: each (measurable) volume is assigned a 
non-negative number, namely the total intensity of radiation scattered into this vol- 
ume. In this article, we will only meet the simple case that v is a Dirac comb u s . 
However, already the convolution with a function of compact support (a "profile" of 
the scatterer) shows why the more general setting is useful. Let us summarize this by 
the following result, which is a combination of |12|, Prop. 2.2] and ]12[ Prop. 3.3]. 

Proposition 7 Let v be a translation bounded measure. If its natural autocorrela- 
tion r y u exists, it is a translation bounded (hence tempered) positive definite measure. 
Furthermore, % is then a positive measure and also translation bounded. 

Decomposition of measures 

A pure point of a measure v is a point x e M. n with u({x}) ^ 0. Since v is a regular 
Borel measure, it has at most countably many pure points and the sum of |y({cc})| 
over the pure points in any compact set K is convergent, and the pure points alone 
serve to define a measure v vv called the pure point part of v. Thus v has a unique 
decomposition as 

v = v p p + Vc ( 44 ) 

where v c := v — v is the so-called continuous partf\ of v, and is characterized by 

having no pure points. A pure point measure is a measure whose continuous part is 0. 

When a tempered measure v is decomposed in this way the components v vv and 

v c are, of course, measures, but not necessarily tempered. For a translation bounded 

6 Note, however, that the definition of the analogous operation to * in |l2f is slightly incorrect in 
that it omits the extra complex conjugation. 

7 The word "continuous" here does not refer to being continuous as a function, since a line measure 
in the plane, for example, is continuous. It refers to the intermediate value property that if there 
are sets A C C with v c (A) < b < v c (C) then there is a set B with A C B C C and v c (B) = b. Some 
authors use the word diffuse instead of continuous and the words atomic or purely discrete instead 
of pure point. 
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measure v, however, the components are both translation bounded and hence tem- 
pered. If v is a positive tempered measure, the decomposition is automatically into 
tempered components. 

Returning now to diffraction patterns, we meet the special situation that %, is a 
•positive measure. Consequently, it decomposes into a pure point part and a continuous 
part, both of which are positive measures. The pure point part is called the Bragg 
spectrum (of the point set S that created it). We say that the diffraction pattern is 
pure point or that S has a pure point diffraction spectrum if its diffraction pattern is 
a pure point measure, i.e. if it is equal to its Bragg spectrum. 

Let us finally mention that, relative to Lebesgue measure, the continuous part u c 
can be further decomposed into an absolutely continuous and a singular continuous 
part, v c = u ac + u sc , see ||, |4j] for details. In our examples, we show that v c vanishes, 
which means that there is neither an absolutely continuous nor a singular continuous 
component present. 

An intuitive derivation of the point spectra 

This section describes a short, intuitive way of calculating the pure point part of the 
diffraction spectrum of the visible points V of a lattice T and also of the /cth-power-free 
numbers F = F^. Its purpose is to give a taste of our later number theoretic methods 
in a simpler setting and also to contrast this intuitively clear calculation with the 
more circuitous route via autocorrelations we take later in rigorously establishing the 
complete diffraction spectra. It seems almost miraculous when the longer method 
eventually reduces to the same simple result. The intuitive method depends on the 
fact that since the autocorrelation 7^ is a volume-normalized convolution of 10 — u v 
with itself, its Fourier transform, the diffraction pattern j u , should be a normalized 
square of Cj. 

In this context, the coefficients of u are usually called amplitudes, even if they 
only exist formally, while those of 7 W are called intensities which relates to the fact 
that they are real and non-negative. The appropriate operation now is to determine 
the intensities as the absolute squares of the corresponding amplitudes of peaks, as 
indicated in the commutative Wiener diagram in Figure |I] where the vertical arrows 
represent the Fourier transform, the upper horizontal arrow the volume-normalized 
convolution, and the lower horizontal arrow taking the absolute squares of amplitudes. 

This observation dates back at least 100 years in optics and is the standard proce- 
dure in diffraction theory [[?]. In the context of diffraction from infinite arrangements, 
it has been made rigorous (at least for the pure point part of the spectrum) by Hof 



FT 




Figure 1: Wiener diagram 



T2| , |13|| . The intuitive approach to the visible points runs through the Wiener diagram 
along the "low road" , via the difficult to interpret Cj. 



Visible points 

We start from 



CO 



v 



Ert 



m )u) 



?nr\{0} 



(45) 



m=l 



which follows from ( |i~5|) since the sum on the right is supported on T \ {0} and the 
amplitude of the peak at x e T \ {0} is 



XvW 



Yl ^( m ) 



(46) 



m|cont(a;) 



Although the sum on the right of fl46|) is finite for each x, these pointwise sums are 
not uniformly convergent, so the sum on the right of (|45"D does not pass the M-test 
of Proposition ||] as a sum of pure point distributions. Nevertheless, it is clear that 
the sum does converge to uj v in the weak*-topology (though not in the local norm 
topology). Taking the Fourier transform of (^) term-by-term we obtain 



UJ, 



Y\ Km) ( —u rt 

m=l x 



/m 



(47) 



where here and in the remainder of this subsection we assume V to have unit density, 
i.e. vol(r) = 1. Since u v is a tempered distribution it has a Fourier transform Cj v 
which is also a tempered distribution, and since the Fourier transform is a continuous 



operator on the space of tempered distributions with the weak*-topology |24| the sum 
on the right of (|47|) converges to Cj v in the weak*-topology (though again it does not 
pass the M-test). This time the sums of the amplitudes at individual points of Qr* 
are infinite sums but uniformly convergent, when n > 2. 
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It is difficult to interpret ( ffTD as a distribution and we suspect that it is not a 
measure. Nevertheless, we can identify the pure points and their amplitudes formally 
(even though they may not be real peaks). 

Define the denominator of a non-zero point x G Qr* by 

den(cc) := gcd{m G Z \ mx G V*} (48) 

(it is the smallest positive integer such that mx G T*). Then the sum of the amplitudes 
in (|47| ) at a point x G Qr*\{0} with denominator d is 



Efi(ld) ji{d) ^-\ nil) n(d) T-r / 1 

(ld) n d n 2-^ l n d n J-l I p n 

i=i v y l=1 pfi v F 



(Z,d)=l 

jJ,{d) T-r/ 1 \~' /i((i) 



n(i~) =^n=^r- m 



When n > 2, the pointwise sum of the pure point parts is thus nonzero at all points 
of Qr* with squarefree denominator but is not absolutely locally summable. It is, 
however, locally square summable. The squaring operation gives, for the diffraction 
spectrum 7 W of V, the distribution with a pure point at each point of Qr* with 
squarefree denominator, the peak at such a point with denominator d having intensity 

(>){J(P"-1) 2 ' (50) 

These coefficients are uniformly absolutely locally summable and in fact are the correct 
intensities of the diffraction pattern. This will be proved below in Theorem |[ They 
were first derived, based on similar arguments, in 0. However, one may well ask why 
this works. How can we justify the squaring operation as the appropriate mechanism 
for obtaining the intensities? The intuition is supported as follows. 

Hof 's results 



Two results of Hof 112, Ol are 



Proposition 8 Let v be a translation bounded measure on W" 1 such that v is also a 
translation bounded measure. Then 

0({x}) = lira -±- [ e- 2 —y u(dy) (51) 

K-oo v n R n J BR{a) 

for every x G M n and the limit exists uniformly in a. 
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Proposition 9 Let v be a translation bounded measure with natural autocorrelation 
7 and suppose that for all x e M. n , 

m x := lim -i- f e^** v{dy) , (52) 

exists uniformly in a. Then, for all x, we have 

7({*}) = \m x \ 2 . (53) 

Taken together, these show that if the Fourier transform v of a translation bounded 
tempered distribution v is also translation bounded (and hence can be decomposed 
as a pure point part and a continuous part) then the pure points of the diffraction 
spectrum 7^ are the same as the pure points of the Fourier transform 0, but their 
intensities are the absolute squares of the amplitudes of v. Proposition |S] alone says 
that this continues to be true even when v is not translation bounded, at least as 
long as the formal expressions for the pure point amplitudes, also called Fourier-Bohr 
coefficients, of v (which may now not represent peaks of v in the accepted sense) 
converge uniformly with respect to translation of physical space. We shall say more 
about the status of the method after giving the analogous intuitive derivation for the 
/cth-power-free integers. 

Mh-power-free integers 

The parallel intuitive calculation for the pure point part of the diffraction spectrum 
of the set F = F k of fcth-power-free integers in R goes like this. We have 



00 



u F -- ^^H^zyo), (54) 



m=l 

hence 



S^H^S/m*- 1 )- (55) 

Let dk be the smallest positive integer such that d\ is divisible by d. We note that 
dk is a divisor of d but is divisible by every prime factor of d and that d \ m h if and 
only if dk \ m. Then the formal sum of the amplitudes of the peaks of (|55|) at a point 
x G Q with denominator d is 

v^ M^fc) M4) ^> MO Mjfc) TT (-, 1 V 1 K<T) TT l /c fi N 

2^{ld k f "^ d\ £j i k " C(k)d k k ^\ P k ) ' ((k)Up k -V { ] 

(l,d)=l 
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since, when dk is squarefree, we have dk = d* = YlpuP (^ ne squarefree kernel of d) 
and d'l = Yi P \dP k - Consequently, the pure point part of the diffraction spectrum of F 
has peaks at the points of Q with (k + l)th-power-free denominator, and its intensity 
at a such a point, with denominator d, is 



( k )li( p k_ 1 )2- ( 57 ) 



This is a pure point distribution and also agrees with the result of Theorem [5] below. 

Status of the method 

In the remaining sections of this paper, which contain the heart of the proofs, we 
completely avoid Hof 's formula d52|) (which in any case could only determine the pure 
point part of the spectrum without any assurance that there is no continuous part). 
What prevents us from making our derivation of the pure point part of the spectrum 
rigorous by citing Proposition ^] is that for our examples the limit in (^) is not 
uniform in a. This is most easily seen by noting that our point sets have arbitrarily 
large holes, so that no matter how large R is there will be an a such that -Br(O) + a 
lies entirely within a hole. However, large holes are very sparse, so the limit is, in 
a sense, "nearly uniform". For x (jL (QT* the limit ( |52"D can be shown to be (as it 
should be) but again is almost certainly not uniform in a. 

The fact that our examples are outside the domain where Hof's result applies 
underscores the fact that we are in new territory not only for calculating the continuous 
part of the diffraction spectrum but even for calculating the diffraction peaks. In other 
words: the answer to the question which distributions of matter diffract is still largely 
unknown. 

Autocorrelation of the visible points 

The remainder of this paper consists of a rigorous derivation of the results of the 
previous section via the "high road" of the autocorrelation. So we begin by using an 
elaboration of the kind of Mobius inversion argument used in proving Proposition |6] 
to calculate the autocorrelation of the visible points V of a lattice V. 

Theorem 1 For n > 3, the natural autocorrelation 7 of the set V of visible points 
of a lattice T exists and is supported on T, the weight of a point a G T in the auto- 
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correlation of V being given by 

[a) = dens(r) £(n) J] (l + -L^\ , (58) 



with error term (as defined just after (Tl)) equal to 0(1/ R), where the implied constant 
depends on a as well as on T. 

PROOF: By Prop. ||, V — V = T, so the autocorrelation of V (if it exists) is supported 
on T. The weight of a point a G T in the autocorrelation of V is the limit as R —* oo 

£f E 1 (59) 

\x\,\x—a\<R 
x,x—adV 

and by Lemma |I] the existence of this limit for each a e V is sufficient to ensure the 
existence of the autocorrelation. 

It is convenient to drop the condition \x — a\ < R in (|59|), which then becomes 



^E 1 ' W 



v n R 

\x\<R 

X,X—a£:V 

The difference between these sums is 0(1/ R), due to the extra lattice points x within 
a constant distance \a\ of the boundary of -Br(O) that are included in the latter. We 
note that the latter sum has a natural geometric interpretation as the proportion 
of integer points in a large ball that are visible both from the origin and from the 
viewpoint a. By (|T5j) and (|S|), ( |B0D can be expressed as 

^ E E moEa*h = ^E E momhE 1 ' ( 61 ) 

£cer\{0,a} i|cont(£c) m|cont(a;-a) 1<KS l<m<S xeF\{0,a} 

\x\<R x=0 (mod IT) 

x=a (mod mr) 
|as|<-R 

where S = (R+\a\) / L(T) . Collecting together terms with the same value old = (l,m), 
noting that all solutions x of the congruences belong to <ir and that the congruences 
have no solution at all unless a e dT, and putting I' = l/d, m! = m/d, x' = x/d, 
a' = a/d, we obtain 

^ru E E E ^I'dMm'd) E 1 • (62) 

d[cont(o) l<l'<S/d l<m'<S/d x'er\{0,a'} 

(l',m')=l x'=0 (mod IT) 

x'=a! (mod mT) 
\x'\<R/d 
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Since I' and m' are bound variables of summation and x 1 and a' will not be referred 
to again, we can drop the dashes: from now on I and m are the new /' and m' but a 
is the original a. 

By Propositions and [l] (with ImT in place of V) the inmost sum is 

& ™ {T ^{S +o {itT +o{i) (63) 

These three terms give a main term and two error terms in (|62|). 
The first error term is majorized by 

JjyJ y l \ fO(l/R), ifn>3, 

V^ i<£* /(n_1) i<t^5 m(n_1 7 " I ^((logS) 2 /^), if n = 2, l ° j 

since the sums are convergent when n — 1 > 2 and increase logarithmically when 
n- 1 = 1. 

The second error term is majorized by 0(S 2 /R n ) = 0(1/ R n ~ 2 ). So when n > 3 
both error terms are 0(1/ R) (since S = O(R)) and tend to as R — > oo. 

The main term is (with D = dens(r)) 

2-^2-^2^ (dlm) n 2^, 2^ 2^ (dlm) n ^ ' 

d|cont(a) l<KS/d l<m<S/d d|cont(a) \<KS/d l<m<S/d 

(l,m)=l (l,d)=l (m,d)=l 

(l,m)=l 

since /i(ld) is ii(l)n(d) when (l,d) = 1 and otherwise, and similarly for fj,(md), 

= dv^v y # (66) 

d|cont(a) l<KS/d \<m<S/d 

(l,d)=l (m,d)=l 

since fi(l)fi(m) is fi(lm) when (Z,m) = 1 and otherwise 

_ p y ^y.MrMO ^^^^ (6?) 

d|cont(a) r=l 

(r,d)=l 

since the double sum is absolutely convergent. The difference between this limit and 
the partial sum (|66|) is 0(1/ i?™" 1 ), so falls within the error term estimate 0(1/ R). 
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Using QT^D, the expression for the limit flST| ) can be rearranged as 

"E £11 (i-£)=»«»)£ £11 (i-£)" < 68 > 



<2|cont(a) P/T d d|cont(a) pM 

squarefree d squarefree 



»««) n ( i + £ - 1) "') = »««) n('v 



, , , , , ,2 

p|cont(a) \ / p|cont(a) 

This completes the proof. □ 

To establish the existence of the autocorrelation function for the visible points 

when n = 2, we need to reduce the second error term in the above proof. This we 



do by modifying the argument slightly, first using the characteristic function (p~5|) to 
replace the constraint on x — a only, then discarding large values of m from the sum 
before using (|15| ) again to replace the constraint on x. 



Theorem 2 Theorem^ holds for n = 2 with the error term increased to 0(1/ yR). 
PROOF: Using ( |T5| ) to replace the constraint on x — a, ( J6QJ ) becomes 

^E E m™) = ^ E m^)E !> ( 69 ) 

jeGV m|cont(a:— a) l<m<S xSV 

Xy^a x^a 

\x\<R x=a (mod mr) 

\x\<R 

where S is as before. The inner sum is trivially 0(R 2 /m 2 ) (since m < S = 0(R)) so 
the contribution to fl60"|) from the terms with m > \^R is 

°{ E ^) =0(l/v^). (70) 



For the terms in fl69|) with m < yR we use the characteristic function (|T5| ) to replace 
the constraint on x, obtaining 

^ E ^E Emo = J^E E mom^E 1 - ( 71 ) 

l<m<-/R ccer\{0,a} l\cont(x) 1<KS l<m<VR ccer\{0,a} 

x=a (mod mr) x=0 (mod W) 

|a;|<_R x=a (mod mr) 
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This is now identical to fl6~T|), except that the second S is replaced by v-R, and as 
before contributes to (|60|) a main term that tends to 

p|cont(a) 



as i? — > oo and two error terms 0((logi?) 2 /i?) and 0(l/y/R). The difference between 
the main term and its limit is also 0(1/ yR). So the total error term is majorized by 
0(1/ VR). D 

Note that this modified argument does not reduce the error term when n > 3, 
which contains a term 0(1/ R) arising from the boundary of the ball when I = m = 1. 
For n = 2, a similar but more complicated argument gives an improved error term 
0(i?~ 3 / 4 (log-R) c ) for some constant c, but we do not need this here. 

Remark Theorems |IJand[^ show that the weight w(a) in the autocorrelation of 
V (also called the autocorrelation coefficient^) depends only on the content of a and 
the density of I\ This enables us to calculate how the application of a non-singular 
linear transformation T to V affects the autocorrelation of V . Clearly T preserves 
content, in the sense that, for any lattice V and seT, the content of Tx as a vector 
of TT is equal to the content of a; as a vector of T. In particular, TV is the set of 



visible points of TT. The autocorrelation of TV is supported on TT and (58) shows 

that 

,_ dens(Tr) 1 

WTviTa) = i^y Wv{a) = Wr\ Wvia) > (73) 

where the suffix on w indicates which autocorrelation it is associated with. Thus the 
autocorrelation of TV is | det T| _1 times the T-image of the autocorrelation of V. 
The direct way of calculating wtv(Tcl) from fl53p is 



<"-< To) = fc^ £ '• (74 > 

n x,x-Ta£RBC\TV 

where B = -Bi(O) is the unit ball in M. n , and substituting this in (|73|) gives 

detTl v-^ , 1 

> 1 = lim — — - 

*-^ R^oovol(£,_ 

y^-aeRT^BnV y ' y,y-a&REnV 

8 Note that w(a) can also be interpreted as the density of lattice points that are simultaneously 
visible both from the origin and from the lattice point a. Extending this to the condition of simul- 
taneous visibility from an arbitrary (but finite) set of points of T results in higher order correlation 
coefficients, see Hi for some recent results. 



Wy(a) = lim -r- 1 > 1 = lim — -r— > 1, (75) 

VK ' R^oo Vn R n *-*> R->covo\(E)R n *-^_ 
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where E = T~ l B is an ellipsoid. Hence averaging over the expanding ellipsoid RE 
gives the same value for the autocorrelation of V as using the natural density and 
averaging over an expanding ball. The choice of ellipsoid here is completely arbitrary, 
since for any E there is a linear transformation T with TE = B. In the Appendix 
this is generalized further and we show how to replace E by any bounded measurable 
region (not necessarily centred at 0) with finite (n — l)-dimensional surface area. 

Diffraction spectrum of the visible points 

The final step in obtaining the diffraction spectrum of the visible points V is to take 
the Fourier transform of the autocorrelation of V. 

Theorem 3 The diffraction spectrum of the set of visible points of an n- dimensional 
lattice T (with n > 2) exists and is a pure point measure which is concentrated on the 
set of points in QT* with squarefree denominator and whose intensity at a point with 
such a denominator q is given by 

dens(r) 2 -,-,- 1 



( 2 (n) -LI (pn 

P\9 



n^rr (76) 



This measure can also be represented as 

oo 1 

densmW^n-^^— K/,> (77) 



p 2n _ 2 p n, 
d=l P\d 



d squarefree 

a weak* -convergent sum of Dirac combs. 



PROOF: Let 7 be the autocorrelation of V. The right hand side of (68[) can be 
expressed in the form (again, with T> = dens(r)) 

^a) = ^H£-n(l--) ■ (78) 

d=l P\d 

d squarefree 
aedT 



So by Theorems |I] and |2| and Lemma 

7 



00 , n\-l 



d=\ P\d 

d squarefree 
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Since ||u; dr || loc = 0(1) and the coefficient of u dr is 0{l/d n ) when n > 2, this 
sum of tempered distributions is convergent in the weak*-topology by Lemma ||] (and 
this is easy to see by other means in this case, where the resulting sum is uniformly 
discrete). Its term-by-term Fourier transform is 

00 1 / 9 \ _1 

^WE^II 1 ^ "TV* (80) 

d=l P\d V V J 

d squarefree 

which weak*-converges to the diffraction spectrum of V, since the Fourier transform 
operator is weak*-continuous. Since ||u; r *; d || loc = 0(d n ) and the coefficient of 0J v */ d is 
0(l/d 2n ), Lemma tells us that the weak*-sum is a translation bounded pure point 
measure equal to the pointwise sum of its terms. 

This establishes the series form (|77| ) for the diffraction spectrum. The explicit 
values of the intensities can now be evaluated quite simply. Let p be a point in Qr* 
with denominator q. We can suppose that q is square-free, since otherwise there is no 
contribution to ( ff7| ) at all. The terms in (|77|) that contribute to the intensity at p are 



those with d = mq {m squarefree in Z + and prime to q), so the intensity at p is 

1 °° 1 

^wn-^^z n^^- ( 8i ) 

P\l m=\ P\m 

m squarefree 
(m,q)=l 



This simplifies to 



^^II^rra^lK 1 



1 



p2n _ 2pn J. A V ' p 2n — 2p n 

p\q pYi 



- 1- „ , 1 s 3 , 2 n - 



^^ ' J. J. T)2n\ rgTl J J-J-V T) n / \ r>" 

p\q p/q 

V 2 t-i- 1 / 1 \-2 

p\1 

(using the Euler products in ([TJ) and (pTQ|) ) which agrees with ([76|) . D 

An explicit example of the diffraction (for T = Z 2 ) is shown in 0, and compared 
with an optical experiment. The diffraction image in this case is both Z^-symmetric 
and GL(2, Z) invariant which results in a beautiful image with a rather unusual sym- 
metry structure, reminiscent of self-similar patterns common in fractals. 
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/cth-power-free numbers 

In this section we derive the diffraction spectrum of the 1-dimensional set F consisting 
of the /cth-power-free numbers in Z. Again, this has arbitrarily long gaps but it 
nevertheless has a pure point diffraction spectrum. The proof of the second assertion 
closely parallels the corresponding proof for the visible points V in the previous two 
sections, with the parameter k for F playing the role of the dimension n of V in the 
formalism. There are some differences of detail, however, particularly with the error 
terms. 

Let us note here that the results on the autocorrelation derived below also follow 
from |22| (k = 2) and from [|19 , Thm. 1], where also more general correlation functions 



have been derived. The error term given in [|1| is 0(R~ 1+2 ^ k+1 ^), which is slightly 
better than the error term we derive in Theorem f|. We include our proof to make this 
paper self-contained and to show up the close parallel between visible lattice points 
and fc-free numbers. 

Proposition 10 F is uniformly discrete, but has gaps of arbitrary length. Moreover, 
for any L > the set of gaps of length at least L has positive density. 

PROOF: The uniform discreteness is trivial since F C Z. Choose L integer moduli 
m 1 , . . . ,m L that are > 1 and coprime in pairs (for example, the first L primes). By 
the Chinese Remainder Theorem there is an integer N with 

N=-j + l (mod m)) for j = 1, . . . , L. (83) 

Now for x = N (mod m k m k . . . m k L ) we have m\ \ x, m\ \ (x + 1),. . . , m k L \ (x + L — 1), 
so none of the numbers x,x+l,...,x + L — l is /cth-power-free. The numbers x have 
density {m l m 2 . . . m L )~ k . □ 

This argument gives a distance of the order L kL between gaps of length L, so again 
long gaps can be expected to be extremely sparse. Nevertheless, gaps of length L have 
a definite frequency (its expression in terms of Dirichlet series can be extracted from 
|JT9||). It is interesting to note that the corresponding distribution is not Poissonian, 



but that these frequencies decline faster than exponentially in L |T5|, 



Proposition 11 For k > 2, the kth-power-free integers F have a natural density 
given by 

dens(F) = -L (84) 

with error term 0{R^ l+<yl / k ' > ). 
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At least for the squarefree numbers, this is again a standard example of Mobius 
inversion, see [14, Thm. 6.6.1]. 

Proof: The natural density of F is the limit as R — > oo of 



2R 

\x\<R 



which by (O) is 



-L E 5>(-) = ^ E "WE 1 (86) 

[x[-CR m fe |a; l<m<K 1 / fc |x|<_R 

xj^O x=0 (mod m k ) 

x^O 

sr aw i n( Rl,k \ 

^ m k V R J' 

l<m<R 1 ' k 

since the inner sum on the right of (|86| ) is 2R/m k + 0(1). This last expression tends 
to l/((k) as R — > oo when n > 2, with the errors due to the tail of the sum and to 
the explicit error term both 0(R~ 1+ ^ k '). □ 

Theorem 4 For k > 2, the natural autocorrelation of the set F of kth-power-free 
integers exists and is supported on Z ; the weight of an integer a in the autocorrelation 
of V being given by 



«wn( i +^). ( 87 > 



w(a) 

p k \a 

with error term 0(R~^ 1 ~^ 1 ' k '' ), where the implied constant depends on a. 

PROOF: Clearly the autocorrelation is supported on Z. The weight of an integer a in 
the autocorrelation of F is the limit as R —>■ oo of 



h E ' «»> 



2R 

\x\,\x—a\<R 
x.x—a fc-free 



and, as in the proof of Theorem [l], to show the existence of the autocorrelation it is 
enough to show that this limit exists for every a. Again as in the proof of Theorem [3], 
we can drop the condition \x — a\ < R from the sum to obtain 



^5 E ' <89) 



2R 

\x\<R 
x.x—a fc-frcc 
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with error 0(1/R). 

This can be evaluated by an argument exactly parallel to that for the visible points, 
using the characteristic function fll6|) in place of ( |l5f) . We use the argument in the 
modified form we used in the proof of Theorem || since, in the context of fcth-power- 
free numbers, it gives an improved error term in all cases owing to the fact that, for a 
1-dimensional set, errors due to the boundary of a large region are trivial. (Using the 
unmodified form would give error term 0(R~ 1+ ( 2 / k ^) instead of the 0{R~ l+<K2 ^ k ^^ l ^ k ^ ) 
we obtain here. As with the visible points, this would not be small enough to establish 
the existence of the limit for the squarefree numbers k = 2 — the paradigm case.) 

Using the characteristic function fll6|) to replace the constraint on x — a, (JS9|) 
becomes 

2R E E "H = 2l E M™)£l, 00) 

\x\<R m k \{x-a) \<m<S 1 l k \x\<R 

x fc-free x fc-free 

x^a x^a 

x=a (mod m k ) 

where S = R + \a\. The inner sum is trivially 0(R/m k ) (since m k < S) so the 
contribution to fl39|) from the terms with m > T, where T is a parameter tending to 
infinity with R to be chosen later, is 

°(E^) =o(r 1 -*). (9i) 

m>T 

For the terms in fl90f) with m<Twe use the characteristic function ( |16|) to replace 
the constraint on x, obtaining 

i E M m )E Emo = ^ E E mom^E 1 - ( 92 ) 

l<m<T |x|<_R i fc |a; IKKR 1 ^ \<m<T \x\<R 

x^0,a x^0,a 

x=a (mod m k ) x=0 (mod l k ) 

x=a (mod m k ) 

Collecting together terms with the same value of d = (l,Tn), noting that all solutions 
x of the congruences are divisible by d k and that the congruences have no solution 
unless d k \ a, and putting V = l/d, m' = m/d, x' = x/d k , a' = a/d k , we obtain 

^E E E m^m^E 1 - ( 93 ) 

d k \a l<V<R 1 / k /d l<m'<T/d \x'\<R/d k 

(l',m')=l x'^0,a' 

x'=0 (mod l' k ) 
x'=a' (mod m' k ) 
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As in the proof of Theorem [TJ, we can drop the dashes from now on. 
By the Chinese Remainder Theorem the inner sum is 



(dim) 

giving a main term and another error term in 
The error term is majorized by 0(R~ l+<yl / k ^T). 
The main term is 

Ev^ v^ ji{l<£)n{md) 

Z_^ Z_^ (dlm) k 

d k \a l<l<R 1 / k /d l<m<T/d 



op 

+ 0(1), (94) 



- E E E "'"?y ( 95 > 

d*|a l<KR l / k /d l<m<T/d 
(l,d)=l (m,d)=l 

(l,rn)=l 

since fi(ld) is n(l)fi(d) when (/,<i) = 1 and otherwise, and similarly for fi(md), 

v^ l* 2 (d) y^ y^ fj,(lm) 

2^i $ 2^i 2-^i (lra) k 

d k \a l<l<R 1 ^ k /d l<m<T/d 

(l,d)=l (m,d)=l 

since fi(l)fi(m) is fi(lm) when (7, m) = 1 and otherwise 

- E^E^ -*-« w 

d k \a r=l 

(r,d)=l 

since the double sum is absolutely convergent. The difference between this limit 
and the partial sum (|96p is 0(1/R~( k ~ l " k ) + 0(T~( fe-1 )) and each of these terms is 
subsumed by one of the error estimates we already have. 
This last expression is 

Ein^HwEinO-^)" 1 «*> 

d k \a Pj« V ^ 7 d fc |a PM 7 

d squarefree d squarefree 

p fc |a \ ^ V ^ 7 / p fc |a V ^ 
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Finally, we choose T = i?( 1//fc ) W k ) 2 ^ making both error terms 0(R 1+ ( 2 / fe ) (V^) 2 ). □ 
We can now fill a small gap in our analogy by the following result. 

Corollary 1 For k > 2, we have F — F = Z. 

Proof: Since F — F C Z, 7^ is clearly supported on Z. From Theorem ||, we get 
w(a) > for all a e Z, hence also Z C F - F. □ 

Let us finally describe the diffraction of /cth-power-free integers. 

Theorem 5 The diffraction spectrum of the set F of kth-power-free integers exists 
and is a pure point measure. It is supported on the set of numbers a/q e Q with q 
being (k+ 1) -power- free. The diffraction intensity at a point with such a denominator 
q is 

727fcyll(p*_i)2- (") 

This measure can also be represented as 

^)E(n^-vK^ (loo) 

d=\ p\d 
d squarefree 

a weak* -convergent sum of Dirac combs. 

PROOF: As in the proof of Theorem |^, it follows from ( p8]) that 



^E^n^-^j ^' ( iQi ) 



x — '■ ■'■ rr / - \ 

7 

d=l P\d 

d squarefree 



where 7 is the autocorrelation of F and this sum of tempered distributions is conver- 
gent in the || • || loc -topology for k > 2. Its term- by-term Fourier transform is 

00 1 / 2 \ _1 

Csll 1 -? "w*> ( 102 ) 

d=l P\d V y J 

d squarefree 

where the sum of the local norms of the terms is convergent, so by Lemma ^| ( |102j ) 
weak*-converges to a translation bounded pure point measure equal to the pointwise 
sum of its terms. Since the Fourier transform operator is weak*-continuous this pure 
point distribution is the diffraction spectrum of F. 
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This establishes ( |10(J| ). To obtain ( p3| ) we note that the terms in ( |100| ) that con 



tribute to the intensity at a/q (where q has to be (k + l)-free) are those with d = mq*, 
where q* is the square-free kernel of q and m G Z + is square-free and prime to q. Thus 
the intensity at a/q is 

1 °° 1 

^n^r^E n p -^^- ( io3 ) 

p\q m=l p\ m 

m squarefree 
(m,g)=l 

This is (|81~D without the factor V 2 and with n replaced by k, so reduces to (|99|) . D 

Further connections and directions 

Above, we have emphasized that the sets of visible lattice points Vp and the set of 



fc-th-power-free numbers F k differ from any regular model set (see [TJ, |(], |27|, |28| for 
definitions and properties) by a set of positive density, suggesting that they cannot be 
obtained from the cut-and-project construction in any natural way. However, there is 
a way of obtaining these sets by cut-and-project using the rational adeles instead of 
Euclidean space as the hyperspace (or embedding space) and using windows which, 
although they have empty interior, are quite natural sets in this context. From this 
point of view, Vr and Fk are "super-singular" model sets. This comes about because 
these sets are the result of sieving over primes: F^, for example, is what remains of Z 
after removing the zero residue class mod p k for each p. 

The cut-and-project construction can be pictured like this. 

B <?- A ^ C 

U (104) 

L 

Usually the hyperspace A is M n+m , with the physical space B = M. n and the internal 
space C = M. m being complementary subspaces having n and 7r int as the associated 
projections with kernels C and B respectively, and L being a lattice in A whose images 
in B and C are dense. A bounded acceptance domain or window Q C C is chosen 
and the model set in B is 

A = A(fi) := {n(x) \xEL, n- mt (x) E Q} . (105) 

Then A is a uniformly discrete set, and also relatively dense if Q has non-empty 
interior. When B and C are orthogonal and Q is a simple region we have 

dens(A) = vol(fi)dens(L) = vol(fi)/vol(.F L ), (106) 
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where T L is a fundamental region of L. This construction admits the generalization 
where A is allowed to be an arbitrary locally compact Abelian group and L a discrete 
subgroup with A/L compact (so that L plays the role of a lattice), see [^] and refer- 



ences therein. It is well known [[18|, ^, ^7]] that, even in this more general situation, 
A(O) is relatively dense if Q has non-empty interior and possesses a uniform density 



if, in addition, <9f2 has Haar measure 0, see [|I3 , |27f for details. 

To obtain the /cth-power-free numbers, we take A to be the ring of rational adeles 
Aq. Let us explain what this is. (For the technical details and background to p-adic 
numbers and adeles see 0, which deals with the case of a general algebraic number 
field in place of Q.) Given a prime p G Z + the p-adic valuation on Q is defined, 
for q ^ 0, by \q\ p = p~ r , where q = p r a/b with r G Z and a, b not divisible by p. 
This satisfies not only the triangle inequality but the stronger inequality \qi + q2\ P < 
max{|gi| p , l^lp}- (Such valuations are called non- Archimedean.) The field Q p of p- 
adic numbers is the completion of Q with respect to | • \ p (analogous to K being the 
completion of Q with respect to the ordinary absolute value) and is locally compact 
in the p-adic topology. In view of the strong triangle inequality, the p-adic numbers a 
with \a\ p < 1 form a ring Z p , called the p-adic integers. (It is the closure of Z in Q p .) 
For any p-adic number b G" Z p the p-adic open ball {a \ \a — b\ p < 1} is disjoint from 
Z p . Hence Z p is closed in Q p and, since it is bounded, also compact. But Z p is the 
disjoint union of the p open balls {a \ \a — i\ p < 1} with i = 1, . . . ,p, hence is also open 
in Q p . As a locally compact Abelian group under addition, Q p has a Haar measure, 
unique up to a multiplicative constant, which can be normalized so that vol(Z p ) = 1. 

Any fixed power of a non- Archimedean valuation is also a valuation, and topolog- 
ically equivalent; but up to this equivalence the p-adic valuations and the ordinary 
absolute value are the only valuations on Q. The rational adele ring Aq is the re- 
stricted direct product with respect to the Z p 's of the completions of Q with respect 
to these valuations, i.e. 

A = la = (a (a )) I a °° G R ' " p G ^ p ' &nd " p G Zp \ (107) 

^ I p for all except finitely many p. J 

and has the restricted product topology for which the sets O^ x n«^p> with O^ 
open in R, O p open in Q p and O p = Z p for all except finitely many p, form a base of 
open sets. As a restricted product of locally compact sets with respect to compact 
sets, Aq is locally compact. 

Now Q embeds in Aq diagonally (i.e. each q G Q can be identified with the adele 
all of whose components are q), and with this identification Q is discrete in Aq and 
Aq/Q is compact. A fundamental region for Q in Aq is [0,1] x Y\ p 1< p , which has 
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volume 1 in the normalized Haar measure on Aq. This gives rise to the interesting 
and natural cut-and-project scheme 

r **- \ ^ ri(^) Qp 

U p (108) 



with componentwise projections and Ylu ) denoting the restricted product. The image 
7r(Q) is, of course, dense in R and the denseness of 7T int (Q) in Yl^p * s equivalent to 
the Strong Approximation Theorem |], §15]. If we choose for the window the closed 
and open set Q := Ylp^p we obtain the model set A(Q) =ZgR. Since vol(fi) = 1 in 
the normalized Haar measure on Yl^p an d dens(Z) = 1 we notice that this satisfies 
( |106| ). If instead we choose f2 := n»(^p \p k ^p) with k > 2 we obtain the "thin" 
model set A(f2) = F^. We say "thin" because, being not relatively dense, Fk is not a 
regular model set. The reason that Schlottmann's result does not apply here is that 
Q has empty interior: Q contains no basic open set because it projects to a proper 
subset of Z p in every non-Archimedean component, whereas basic open sets project 
to Z p in all except finitely many components. Nevertheless, Q has a positive volume 
in the normalized Haar measure on Y[ Qp given by 

vol(O) = l[vo\(Z p \p%) = J](l-^) = _L (109) 

which, in view of Proposition [11], agrees with ( |106| ). This is astonishing and we have 
no explanation for it at present, since merely by translating Q we can not only reduce 
the density of A to but can cause A to vanish altogether. Let i : P — > Z, where 
P is the set of positive primes, be any one-one correspondence, and for each p e P 
choose a p G Z p with a p = i{p) (modp fc ). Then the window Q + a C \\7L V , where 
a = (a p ) G Yl Qp; leads to the empty set under the cut-and-project construction 
because each i G Z is excluded mod p k for the corresponding p. 

The visible points Vr of a lattice T C M n , for n > 2, can be obtained by a similar 
construction, using (Aq)", which is topologically isomorphic to the restricted product 
with respect to Z" 

A := K-xJ]^) Qp- (HO) 

p 

Choose a basis {6i, . . . , b n } of V and embed Qr in A by 

X = q 1 x 1 + ■■■ + q n x n h-> (x, ((q v . . . , q n ) p )) , (111) 
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where (gi, . . . , q n ) p = (qi, . . . ,q n ) for each p. (In the terminology of commutative 
algebra, this establishes that A is isomorphic to T ®i A Q , the AQ-ification of the 
Z-module T, see Jl^j under the index entry "ification"). Then the image of Qr is a 



discrete subgroup of A and A/QT is compact with T = T v x \\ Z™ as a fundamental 
region, where jF r is a fundamental region of r in M. n . The volume of this fundamental 
region in the normalized Haar measure on A is l/dens(r). This gives the cut-and- 
project scheme 

U v (112) 



with the images of Qr under n and 7r int being dense. If we choose for the window the 
closed and open set Q := Ylp^p we obtain the model set A(fi) =Tcl" which has 
vol(fi) = 1 and satisfies OTUop . If we choose ft := H p {Z™\pZ™) we obtain A(fi) = Vr- 
Again, f2 has empty interior but positive volume, given by 

voK«) = n(i-^) = cM (113) 

so that ( |106|) holds by Proposition |6|. Again, translating Q can cause the model set 
to vanish. A translation a. 6 FJ Q™ that does this can be constructed as follows: let 
i : P — > Z n be any one-one correspondence and for each p choose cx. p = i(p) (mod p); 
then take a = (a p ). 

These adelic constructions are nothing more than number theoretic sieves which 
exclude certain residue classes modulo a power of each prime. The choice of window 
determines the residue classes to be retained and translating the window changes the 
set of residue classes without changing their number. Clearly, many other examples 
of the same type can be constructed. 

If Vp and F k were regular model sets, we would be able to use standard results 



on model sets to derive their diffractiveness [pq] , but the constructions just described 
are outside the range of current diffraction results. To the best of our knowledge, the 
diffraction of adelic model sets has not been studied, which is why we have had to 
derive our diffraction results from scratch. 

Obviously, one further question is how relevant the examples are. Is the set of 
visible lattice points an isolated example, or more of a paradigm case of a family 
of examples with pure point diffraction? We believe the latter is true, and there is 
one immediate class of examples that springs to mind. If we view the visible points 
of Z 2 as the orbit of (1,0) under the group SX(2,Z), it is clear that this orbit, in 
general, will split into several pieces when SL(2, Z) is replaced by one of the standard 

37 



congruence subgroups. However, our methods are robust under further congruence 
constraints, and so it is clear that there is a large class of congruence subsets of the 
visible points which will have interesting diffraction spectra. 
We hope to report on some examples soon. 

Summary 

We have demonstrated that the set of visible lattice points in dimensions n > 2 and 
the set of /cth-power-free integers with k > 2 both possess well-defined autocorrela- 
tions and pure point diffraction spectra. We have shown how to replace the intuitive 
but incomplete derivation of the spectra by a rigorous number theoretical argument 
based on explicit computation of the autocorrelation followed by Fourier inversion 
and repeated application of the Poisson summation formula. 
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Appendix 



This paper has required several calculations of densities of discrete point sets, for 
which we have used natural density. The natural density is in fact nowhere near as 
natural as one would expect from its name. Here we briefly examine the drawbacks 
of the natural density and also the uniform density, and suggest a simply stated 
definition of density, intermediate between natural density and uniform density, that 
is applicable to sets of the kind we encounter here. A fuller exposition will appear 
elsewhere 



As an aid to describing the behaviour of densities we begin by introducing some 
standard measurements of measurable sets in R n . 

Measurements of measurable sets 

For a bounded, measurable set 1Z in IR n we denote by V(JZ) the volume of 1Z (that 
is, its n-dimensional measure). Given e > we denote by dJZ the set of points in 
W 1 whose distance from the boundary of 1Z is less than e. Then dJZ is open, hence 
measurable, so also has a volume. We temporarily call a set 1Z good if 

(a) 1Z is bounded and measurable, and 

(b) V(d e 1Z)/e is bounded as e — ► 0. 
For good sets we can define the functions 

(114) 
(115) 

(116) 

(117) 

The following relations are immediate from these definitions: 

R < M, V < v n R n } RS > 2 n ' 1 V, S > v n R n - 1 /2, (118) 

where v n is the volume of the unit ball in R n . 

Each of these functions is invariant under Euclidean transformations and scales 
homogeneously with dilation of 1Z: M and R proportionally to the dilation factor A, 
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V = 


= V(1Z), 




M = 


= M(K) = 


- sup \x\, 
xen 


R = 


= R(K) = 


inf sup x — c , 

ceiR n x&n 


S = 


- S(K) = 


sup {V(d e 7Z)/2e}. 

0<e<R 



V proportionally to A n and S proportionally to A n_1 . They are affected by a non- 
singular affine transformation A = T + t (where T is linear and t is a translation) as 
follows: 



V{A{K)) = \ det T\V(K), 




(119) 


M{A{1Z)) < \\T\\M(R) + \t\, 




(120) 


R(A(K)) < \\T\\R(1Z), 




(121) 


S(A(1Z)) < 2 n |detT| UTII"- 1 !^- 1 ! 


\ n S(lZ). 


(122) 



Of these functions, V, M and R are straightforward, representing the volume, 
the maximum distance from the origin and the circumradius of 7Z, respectively. The 
function S is a substitute for the surface area of 1Z and in fact is always greater than 
or equal to the (n — l)-dimensional measure of the boundary (97?.. The condition e < R 
in its definition is arbitrary, but ensures the scaling of S under dilations mentioned 
above. (Using any fixed multiple of R would have the same effect.) 

These functions extend to arbitrary measurable sets if we define M = R = oo 
when 1Z is unbounded (in which case V may or may not be infinite) and S = oo when 
(b) fails. We note that V(d € lZ) is infinite for every e when 1Z is an unbounded set of 
finite volume, so S = oo for every set of finite volume that is not good. 

Densities of discrete sets 

Let X be a locally finite set of points in R™, that is, every bounded region of R" 
contains only finitely many points of X. Here we compare ways of defining a density 
for such a set X. 

Uniform density 

X has uniform density T> = T>(X) if 

\{xeX | \x-c\ <R}\ = Vv n R n + o(R n ) asi?^oo (123) 

uniformly in c. 

It is immediate from this definition that 

(i) uniform density is invariant under Euclidean transformations. 
A less immediate consequence is that if 1Z is an arbitrary measurable set in R n and 
f{V) is any function of V (no matter how slowly increasing) that tends to infinity 
with V then 

\1ZnX\ = VV(1Z) + o(v(K)\ +0(f{V{K))S{K)\ as V(K) -> oo. (124) 
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A consequence of ( |124j ) is 
(ii) uniform density is independent of shape; 

that is, if X has uniform density T> then ( |123| ) continues to hold, with the same value 
of T>, when the sphere \x — c\ < R is replaced by an expanding set of any other shape 
and v n is replaced by the volume of the set in the family that has R = 1 . 

Another consequence of (|124|) (using ( |119| ) and (|122|) to control the error terms) 
is 

(Hi) uniform density varies as the reciprocal of the determinant under non-singular 
affine transformations. 

All crystals and most quasicrystals possess a uniform density. However, the set 
V of the visible points of a lattice V does not. In fact Proposition 6 shows that the 
number of points of V in a ball of radius R with centre is (dens(T) / ((n))v n R n +o(R n ) , 
whereas Proposition 5 shows that there are arbitrarily large balls that contain no 
points of V. Similarly, the set Fj, of /cth-power-free integers in R does not possess a 
uniform density either. Consequently uniform density is of no use for the questions 
considered in this paper — in fact none of the sets whose densities we require possesses 
a uniform density. 



A definition of density known as van Hove density is implicit in Ch. 2 of PS[, if the 



identical formulae (3.3) and (3.8) of |£6|, Ch. 2] are regarded as defining the density 
of a general potential operator $. This is equivalent to uniform density. The fact 
that uniform density implies van Hove density follows from ( |124f) and the converse 
implication results from choosing the van Hove sets to be balls. 

Natural density 

As defined in (10) for uniformly discrete sets, X has natural density T> = T>(X) if 

\{xeX | \x\ < R}\ = Vv n R n + o(R n ) zsR^oo. (125) 

This is often called asymptotic density in the context of number theory. 

Like uniform density, natural density is 
(i) invariant under Euclidean transformations. 

(Tranlation invariance is easy to derive, though not quite as immediate as for uniform 
density.) However, natural density is not independent of shape and does not transform 
naturally under linear transformations. For example, let X = {(x,y) \ x, y G Z, \y\ < 
2\x\} be the set of points in 1? that lie in the double wedge with angle 2 tan -1 2. The 
natural density of X is (2/n) tan -1 2 = 0.7048 . . . , but the proportion of the integer 
points in the square {(x,y) \ \x\, \y\ < R} that belong to X is 3/4 and the proportion 
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of the integer points in the square {(x, y) \ \x\ + \y\ < R} that belong to X is 2/3. 
So for this set, counting points according to the L 2 , L°° or L 1 norm in R 2 gives three 
distinct answers. The linear map T(x, y) = (2x, y) transforms X into the set of integer 
points with ^-coordinate even in the double wedge {(x, y) \ x,y e Z, \y\ < \x\} with 
angle it/2. This set has density 1/4, which is not a half of (2/7r) tan -1 2, showing that 
natural density does not transform naturally under linear maps. 

For the sets we deal with in this paper natural density does not exhibit such 
pathological behaviour. For example, Proposition 6, Theorem 1 and Theorem 2, which 
establish the densities of certain subsets of a lattice T in M. n , all give a factor dens(r) 
as the only dependency of the density on the lattice, confirming that the density 
varies as the reciprocal of the determinant under linear transformations. However, 
rather than having to verify this in each individual case, as in effect we have done 
here, it would be preferable to have a definition of density with strength intermediate 
between uniform density and natural density for which non-pathological behaviour is 
guaranteed. 

Tied density 

It is possible to define an alternative form of density, intermediate between uniform 
density and natural density, that overcomes these difficulties. This depends on the 
idea of making uniform density less free of the origin by replacing R by M in the error 
term of ( p3|) . We say that X has tied density V = V(X) if for all ceW 1 

\{xeX | \x-c\ <R}\ = Vv n R n + o(M n ) asR^oo, (126) 

where M = R + \c\. 

Again it is easy to see that 
(i) tied density is invariant under Euclidean transformations. 

In a similar way to ( |124j ) it can also be shown that if X has tied density V and TZ is 
an arbitrary measurable set in IR n then 

\nnx\ = vv(n) + o(M(n) n ) +o(m(k)S(k)\ a,sV(n)^oo. (127) 

This is sufficient to establish that 

(ii) tied density is independent of shape, 

and, with the help of ([lT9|) and ( p0|) , that 

(Hi) tied density varies as the reciprocal of the determinant under non-singular affine 

transformations. 



42 



The error term o(M n ) in the definition of tied density has the effect that, while 
balls centred away from the origin are not excluded from consideration, those distant 
from the origin by more than a few radii have little influence. It can be regarded as 
a reminder to experimentalists that, as the radius of a probe is increased, the next 
position of the probe should not be more than a few radii from the previous position. 
There is, nevertheless, a great deal of arbitrariness in the definition of tied density. 



The core idea is to make the error term in ( |12q ) depend on c as well as on R, but any 
kind of dependency on c (provided that for every R the error term tends to infinity 
with c) would have a similar effect. The faster the rate of increase of the error term 
with c the more widely applicable the density definition becomes. In an effort to 
conceal this arbitrariness we have chosen a very simple dependency on c. 

It is clear from these definitions that the existence of uniform density for a set 
X implies the existence of tied density for X, which in turn implies the existence of 
natural density for X, and that each of these densities has a unique value when it is 
defined. 

To establish the existence of tied density for the sets studied in this paper we 
would need to count points x with \x — c\ < R instead of |cc| < R. As a result, 
though the x and y variables in sums still have range R (but no longer centred at the 
origin), the / and m variables have range R + \c\ = M. The effect is that main terms 
remain the same but some i?'s are replaced by M's in error terms. The new error 
terms for the results of this paper when the ball of radius R is centred away from the 
origin are as follows: 

Proposition 6 (n > 3) 0(1/R) + 0(M/R n ), 

Proposition 6 (n = 2) O (log M/R) + 0(M/R 2 ), 

Theorem 1 (n > 3) 0(1/ R) + 0(M 2 /R n ), 

Theorem 2 (n = 2) 0(^JJd/R) + 0(M/R 2 ), 

Proposition 11 (n = 1) 0(M 1 / k /R), 

Theorem 4 (n = 1) 0{R l l k M 1 > k /RM 1 ^ 2 ) + 0{M l / k /R). 

Each of these error terms, when multiplied by R n , is well within the error estimate 
o(M n ) required by the definition (|126|) of tied density. Consequently we could have 
worked with tied densities instead of natural densities throughout. One concrete 
advantage of such an approach would have been that we could then have worked with 
the visible points of cubic lattices only, then used property (iii) of tied density to 
transfer the results to other lattices. 
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